From Cooperative Paramagnet to Neel Order in Y2RU2O7 
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Spin correlations in the pyrochlore antiferromagnet Y2RU2O7 with Curie-Weiss temperature 
Qcw ~ —1100 K and critical temperature Tjv — 77 K were examined through neutron scatter- 
ing. For Tjv < T < Qcw/3 the data show spin relaxation with a rate HF = 1.17(9)fcsT. For 
T < Tjv spectral weight moves to higher energies with substantial changes up to 4 x fcsTjv. For 
T « Tm there is a A = 11(1) meV energy gap and a pronounced spectral maximum at 19.7 meV. 
Throughout the temperature range examined the wave vector dependence of inelastic scattering 
exhibits a broad peak for Qd ~ 3.8 (d is Ru-Ru spacing) consistent with dipolar spin correlations. 
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In the pyrochlore structure, with general formula 
A2B2O7, both the A (trivalent rare earth ions) and B 
(tetravalent transition metal ions) site ions form a lat- 
tice of corner sharing tetrahedra.— Weakly connected and 
non-bipartite, this lattice produces anomalous frustrated 
magnetism when populated by spins with near neighbor 
antiferromagnetic (AFM) interactions.— For classical spin 
vectors the lowest energy states form a manifold char- 
acterized by zero magnetization on every tetrahedroni^ 
Rather than Neel order for T « |0cvf|i pyrochlore mag- 
nets remain disordered to much lower T where the col- 
lective properties are governed by quantum fluctuations, 
thermal fluctuations, and sub-leading intcractions4i^i^ 
Examples of low T states that have been reported are 
spin glasses^, spin liquids^ii^, and magneto-elastically in- 
duced Neel order fiSiii 

Y2RU2O7 is a Mott insulator with 5=1 Ru''+ spin 
and a large negative Curie- Weiss temperature Qcw ~ - 
1100 Kfi^iHiiiii^ Specific heat and susceptibility indicate 
a magnetic phase transition at T/v^ 77 KM^^^^ Diffrac- 
tion measurements show that the transition is to a long 
range ordered (LRO) AFM q ^0 structureJ^ The total 
spin per tetrahedron vanishes so this state is indeed a 
specific member of the low energy AFM pyrochlore man- 
ifold. Here we report a neutron scattering experiment 
probing the anomalous spin dynamics of Y2RU2O7. The 
measurements reveal that AFM spin correlations are es- 
tablished at temperatures at least as high as 300 K and 
evolve little upon cooling. The spin relaxation rate how- 
ever decreases linearly in T according to hT = l.llkBT 
consistent with proximity to quantum criticality. How- 
ever Neel order yields an 11(1) meV gap in the excitation 
spectrum and moves spectral weight to a 19.7 meV peak. 

A powder sample of Y2RU2O7 was synthesised using 
the solid state reaction method. A mixture of Y2O3 and 
Ru02 in proper molar ratio was pre-reacted at 850 °C for 
15 h in air and then ground and pressed to pellets. The 
pellets were sintered at 1000-1200 °C in air with interme- 



diate grindings to improve homogeneity. Powder neutron 
diffraction data were collected using the BTl diffractome- 
ter at NIST. Rietveld analysis confirmed the cubic py- 
rochlore structure above and below Tn with spaccgroup 
FdSni (227) and lattice parameter a = 10.12824(2) A at 
T = 10 K. The analysis also revealed 0.044(2) weight 
% unreacted Y2O3. Inelastic neutron scattering experi- 
ments were performed on the MARI spectrometer at the 
ISIS Facility and the DCS spectrometer at NIST. On 
MARI we used an incident energy E,;= 150 meV for an 
elastic energy resolution of 3.5 mcV. The phonon contri- 
bution to neutron scattering was subtracted by scaling 
the spectrum measured at high wave vectors, where the 
magnetic response is negligible, using the Discus pack- 
age.— This program uses a Monte Carlo method to cal- 
culate a wave vector dependent scaling factor, which 
showed no significant energy dependence. For the DCS 
experiment we used 13.09 meV neutrons, which produced 
an elastic energy resolution of 0.77 meV. For both ex- 
periments count rates were normalised using Y2RU2O7 
Bragg peaks. This produced absolute measurements of 
i{Q,huj) = (7ro)2|fF(0)p25(g,?iw) to an overall scale 
accuracy of 20 %. Here 7= -1.913 and g ~ 2 are gyro- 
magnetic ratios for the neutron and Ru'*"'" respectively, 
ro= 2.82 fm is the classical electron radius, and F{Q) 
is the magnetic form factor. Lacking specific informa- 
tion for Ru^+, we used the Ru"*" form factor for model 
comparison;^ S{Q, Tiio) is the spherically averaged dy- 
namical correlation function^Si 



The magnetic neutron scattering cross section for 
Y2RU2O7 powder is shown in Fig. [1] for T~ 15 K, 50 
K and 90 K. For T > T/v = 77 K the spectra indicate 
short range AFM correlations and a spin relaxation re- 
sponse with a relaxation rate F ksT/h. Upon cooling 
below T/v, spectral weight moves to higher energies, with 
modifications up to 4 x ksT^. Fig.[2]shows that the wave 
vector dependence of the inelastic magnetic scattering is 
all but unaffected by the comprehensive rearrangement 
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FIG. 1: Normalized magnetic neutron scattering intensity, 
I{Q, hw), for a Y2RU2O7 powder sample at T = 90 K (top), 
50 K (middle) and 15 K (bottom). Phonon scattering was 
subtracted as described in the text. The elastic intensity is 
dominated by nuclear scattering. 



of the excitation spectrum. Indeed the Q-dependence of 
scattering changes little below room temperature. These 
data are consistent with frustrated interactions with an 
energy scale set by 10014-' I = 1100 K that establish a low 
energy manifold of near degenerate states with a charac- 
teristic inter-atomic "form factor."— 

The fourier transform of six coUinear antiferromagnet- 
ically correlated spins on the vertexes of a hexagon pro- 
vide a reasonable account of the Q— dependence of in- 
elastic scattering in ZnCr204 over a wide temperature 
range . ^^'^'^ We have calculated a spherical average of this 
function for comparison with the powder data in Fig. [2Ji. 
The close correspondence of calculated peak width and 
position with the data suggest that antiferromagnetically 
correlated hexagons are also prevalent in Y2RU2O7. 

A rigorous understanding of spin correlations in the 
pyrochlore antiferromagnet is now emerging from theo- 
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FIG. 2: Q integrated (a-d) and tiu; integrated (e-h) inelastic 
neutron scattering intensity of Y2RU2O7 taken at T = 200 
K, 90 K, 50 K and 15 K. Ranges of integration are indicated 
in each column. The lines show fits to the data described 
in the text. Lines in (a) and (b) account for the finite sam- 
pling of Q — w space and partially resolved elastic scattering. 
Blue squares in (f ) are normalized energy integrated data from 
DCS covering the energy range from 2 meV to 8 meV. 



retical work. For classical spins {S 00), any spin con- 
figuration on the pyrochlore lattice where all tetrahedra 
have zero net moment is a ground state. Higher order 
spin wave theory {l/S expansion) and the Schwinger- 
boson large N expansion show that quantum fluctuations 
favor states that maximize the number of flippable spins. 
Henley argues that this leads to selection of coUinear spin 
configurations so the low energy physics can be repre- 
sented by an effective Ising model with a collective spin 
axis4 The selected spin configurations further are charac- 
terized by a zero divergence condition for a polarization 
field defined on the dual diamond lattice of bonds. The 
effective Ising model favors AFM hexagons?^ which may 
account for their prevalence in ZnCr204 and Y2RU2O7. 
The corresponding spin correlation function iS(Q) is how- 
ever not the fourier transform of a finite size object as 
it has infinitely sharp features (pinch points) associated 
with the divergence free conditional^ Isakov et al derive 
an analytic expression for iS(Q) from the large N ap- 
proximation, which accurately reproduces results from 
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FIG. 3; T dependence of the spin relaxation rate derived from 
the global fitting analysis described in the text. The solid line 
shows the best straight line fit, TiT = CksT with C = 1.17(9). 



Monte Carlo simulations of the classical spin Heissenberg 
modeli^ The black solid lines in Fig.[2]is the spherical av- 
erage of this expression multiplied by the squared Ru+ 
form factor.— The discrepancy at high Q may result from 
residual phonon contributions to the experimental data 
and/or the ruthenium magnetic form factor employed in 
the model calculation. With the overall amplitude as 
the sole adjustable parameter the large N result provides 
an excellent account of the data. Note though that the 
spherical averaging inherent to powder neutron scatter- 
ing precludes insight on finer details such as the predicted 
pinch points. At elevated temperatures defects should 
appear in the form of monopoles^ that give rise to a fi- 
nite divergence of the polarization field and broadening of 
pinch points^^ which again may not be detectable in pow- 
der data. While the theory is not relevant for the ordered 
state, the close correspondence with the measurements 
there indicate similarities in local dynamic correlations 
across the phase transition. 

Although the Q— dependence of the data does not 
change appreciably below room temperatures Fig. [2] 
shows substantial changes in the fluctuation spectrum. 
Above Tjy the spectra are well described (solid line 
in Fig. EJa-b)) by the following response function: 
x"(lj) = xo^^r/ (r^ + cj^), which is related to the neutron 
scattering spectrum through the fluctuation dissipation 
theorem.— Factorization of the Q and u; dependence of 
the cross section was assumed so F characterizes the lo- 
cal response. The T— dependence of the relaxation rate 
(Fig. [3 is described by hT = CksT where C= 1.17(9). 
Linear in T dependence of F was previously found in sim- 
ulation of classical Heisenberg spins with AFM interac- 
tions on a pyrochlore latticci^ The observation that fcsT 
is the relevant energy scale indicates a spin system close 
to quantum criticality and is consistent with the presence 
of a long length scale from a divergence free polarization 
field. 

For T < Tn a gap appears in the cxcitaton spectra as 




FIG. 4: (a) T dependence of the integrated intensity of the 
(111) Bragg peak (•) which has nuclear and magnetic contri- 
butions, (o) Indicates the postion of the inelastic peak in the 
spectrum, derived by fits described in the text, (b) Results 
of the first moment analysis described in the text, (o) Indi- 
cates the change in spin exchange energy, while the solid line 
denotes the overall change in energy obtained by integrating 
specific heat data.— 

spectral weight moves to higher energy (Figs. [D^c) and 
m^d)). At T = 15 K the neutron data indicate an energy 
gap A = 11(1) meV. This is significantly less than the 
23 meV gap inferred from specific heat data which may 
be dominated by activation of the modes associated with 
the 19.7 meV peakJi A gapped spectrum was confirmed 
by the normalized higher resolution DCS data which in- 
tegrated over energies from 2 mcV to 8 meV show no 
evidence of spin waves emanating from the AFM Bragg 
peaks (solid blue points in Fig. ^h)) as in ZnCr204 
When a phase transition breaks a continuous symmetry, 
Goldstone's theorem assures gapless excitations. The ap- 
parent spin gap in the ordered phase thus indicates axial 
single ion and/or exchange anisotropy. While the gap 
A w 1.7 X UbTn surely is an important aspect of the 
phase transition A is only k, 12% of the ksQcw so spin 
space anisotropy may be less relevant in the high T phase. 

Having examined both the high T critical state and the 
low T ordered state we focus on the transition between 
these phases. Fig. S] compares the T dependence of the 
inelastic peak position to the integrated intensity of the 
(111) Bragg reflection. Apart from probing the staggered 
magnetization squared, the (111) reflection has a small 
nuclear component and therefore does not go to zero for 
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T > T/vji^ The inelastic peak position was obtained by 
fitting a damped harmonic oscillator response^^ to the 
0.8 < Q < 1.4 integrated part of the excitation 
spectra. While the inelastic peak moves to lower energies 
upon heating (see Fig. 4) it is important to note that a 
finite energy peak persists for T = Tn- This indicates 
that short range dynamic correlations at Tn are sufficient 
to maintain a pseudo-gap in the excitation spectrum. 

The shift of spectral weight from a relaxation response 
towards higher energies implies that the energy of the 
spin system, (Hs), is reduced through the phase transi- 
tion. This can be quantified using normalized inelastic 
magnetic neutron scattering data and the first moment 
sum- rule — i21 



3 Lj{l - e-^''^)AS{Q, nLu)du 

2 1 - sinQd/Qd 



(1) 



Here we consider only the isotropic nearest neighbor in- 
teraction between spins separated by a distance d. A 
more complete expression that includes anisotropy terms 
can readily be derivedS^ but the quality of the present 
data do not warrant the added complexity. Limiting the 
integral to 5 to 30 meV and averaging data for 0.8 < 
Q < 1.4 the temperature variation of the spin en- 
ergy A{Hs) relative to the value at T = 90 K is shown as 
open circles in Fig.|3|D. The data can be compared to the 
overall change in energy obtained from the magnetic spe- 
cific heaOi, AEiT) = Gra{r)dT' - j^^^ Cm{T')dr , 
shown as a solid line in Fig.djD. While most of the varia- 
tion in the thermal energy is accounted for by A{TCs) and 
no structural changes have been observed at the transi- 
tion, the data do leave room for an increase in the lattice 
energy below T/v by up to 1.5(5) meV/f.u. Note also that 



spin-phonon coupling was detected through an anomaly 
in the temperature dependence of a Raman active mode 
frequency at Tj 
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While all known pyrochlore antiferromagnets show 
some form of (intrinsic or extrinsic) static correlations for 
sufficiently low T, only Y2RU2O7 and its derivatives fea- 
ture a second order phase transitions to long range Ncel 
order. The apparent lack of involvement of lattice degrees 
of freedom indicates that this is a true magnetic insta- 
bility the understanding of which may require a different 
frame work than for Z11CT2O4. The apparent gap in the 
excitation spectrum indicates that spin space anisotropy 
may be significant. While we detect the magnitude of 
the gap there is no direct information about the magni- 
tude or nature of any spin space anisotropic terms in the 
hamiltonian. The fact that the gap mode is relatively 
sharp in energy for a powder sample indicates that the 
corresponding excitations are weakly dispersive and this 
points to local soft modes pushed to a finite energy on ac- 
count of quantum fluctuations. Indeed the first moment 
analysis oi S{Q, uj) shows that the energy gain associated 
with ordering is intimately related to the development of 
the 19.7 meV peak. The principal challenge posed by our 
experiment is thus to understand the interplay between 
quantum fluctuations and anisotropy terms in establish- 
ing the Neel state spin gap and spin resonance in the 
frustrated pyrochlore AFM. 
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